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We provide a unifying framework for exact, probabilistic, and approximate conversions by local 
operations and classical communication (LOCC) between bipartite states. This framework allows 
us to formulate necessary and sufficient conditions for LOCC conversions from pure states to mixed 
states and it provides necessary conditions for LOCC conversions between mixed states. The central 
idea is the introduction of convex sets for exact, probabilistic, and approximate conversions, which 
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I. INTRODUCTION 



Bipartite entanglement has been shown to be a valuable resource in quantum information theory, allowing two 
^ , parties, Alice and Bob to teleport unknown quantum states [Q, to remotely prepare quantum states using 
classical communication only or to carry out other quantum information processing tasks. It is therefore natural that 
a theory of bipartite entanglement has been developed which attempts to quantify this resource. In this theory, a 
central role is played by the set of quantum operations which are constructed from local quantum operations by the 
parties Alice and Bob and classical communication between Alice and Bob, since these operations cannot enhance the 
quantum correlations in a bipartite state. Here we will call these quantum operations LOCC operations, for Local 
Operations and Classical Communication. The goal of the theory of bipartite entanglement is to obtain a classification 
of states with respect to the set of LOCC operations. 
^ ^' Since the notion of majorization Q is to play a central role in what follows, let us recall some helpful notation. Let 
A and jl be two n-dimensional vectors with real coefficients in decreasing order, A^ > A^ > . . . > A", A* = 1 and 

5Ih ' likewise for jl. We write 

■ 

^1 fc fc 

O^. A^Ai^Vfc = l,...,n^A*>^/i\ (1) 



i=l i=l 



When for all I we have 

n n 
i—l i—l 

we say that pX fl, or in words, fl is weakly supermajorized by pX. For a bipartite pure state {ip) we write A^ 
for the vector of Schmidt coefhcients of {ip) in decreasing order. (If lijj) is unnormalized then A^ is associated with 
the eigenvalues of (TtaIiPXiPI) / {'>P\'>P) ^ again in decreasing order.) In this paper we will call |(V'i|V'2)| and F{p,p') — 
Tr(-\/ pi/^p'pi/^) the square-root-fidelity between states ji/'i) and 1-02) and p and p' respectively. 

Let us begin our study by recapturing what is known for LOCC conversions between pure states. Let \^-!p) be a 
state that we start with, which is characterized by a Schmidt vector fl. 

1. Nielsen A pure state can be obtained by LOCC from the state \ij)p) exactly if and only if 

A^ >■ fl. (3) 

2. Vidal §: A pure state can be obtained by LOCC from the state with probability p if and only if 

pX^ M. (4) 
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The protocol that achieves the optimal probability can be viewed as a two-step protocol 

IV'm) ^exact 10 ^Ma IV'), (5) 

where exact and AAa denote an exact conversion and a unilateral measurement by Alice respectively. The state 
10 is an intermediate state which can be determined by a procedure given in Ref. 

3. Vidal-Jonathan-Nielsen [Q: The optimal / = with which we approximate a state \'ip) by a state j?/;') 

which we obtain by LOCC from is equal to KV'IOI where |0 is the intermediate state in the optimal 
probabilistic conversion of \ipp) into |^). (Note that in cases where the optimal probability is zero, the protocol 
still allows for the construction of the state |0-) 

Our goal in this paper is to extend these results to the domain of mixed states. More specifically, we are interested 
in the 'formation' problem for bipartite mixed states: Given a pure state \'4'p), when can we obtain a density matrix 
p from by LOCC, exactly or with probability pi Similarly, what is the maximum fidelity we can achieve when 
approximating pi The problem of entanglement distillation B for a mixed states is not covered by this formalism. 
What we will find is that the 'formation' problem of mixed states can be completely translated to what is known for 
pure states, i.e. the 3 items listed above. Our results are expressed in four theorems that we present in this section 
and that we will prove later on in the paper. 

Let us first define a restricted class of LOCC operations. We will call a quantum operation an LOCC^ operation 
when it is given by a local measurement /superoperator by Alice consisting of a set of commuting Kraus operators, 
followed by 1-way classical communication from Alice to Bob, and additional local unitary transformations by Alice 
and Bob, which may depend on Alice's measurement outcomes. Similarly, the class of LOCC^ maps is the set of 
LOCC^ maps with Alice and Bob interchanged. 

For every pure state jV'/j) we define three kinds of sets, which we call Si, Si ^ and Si j:, where the first set relates 
to exact conversions, the second to probabilistic conversions and the last one to approximate conversions. These sets 
are defined in the following way 

Definition 1 (Exact) The set is the convex closure of the set of pure states \ip) such that 

Xip >- p. (6) 

The set Si is defined as 

S+ = {£(|V)(V'|) I >C e LOCCf , IV-) G Sp}. (7) 



Note that the sets Sp are generalizations of the Schmidt number sets Sk which were introduced in Ref. A 
Schmidt number set Sk is a set Sp where p are the Schmidt coefficients of a maximally entangled state of Schmidt 
rank k. 

Definition 2 (Probabilistic) The set Sp^p is the convex closure of the set of pure states {ip) such that 

pXi, p. (8) 

The set Si ^ is defined as 

St ^ WWi^l) I ^ e LOCC^ 1^) e Sp,p}. (9) 



In Section II C 2 we will discuss the relation between the sets Si^ and Sk 



Definition 3 (Approximate) The set Spj is the convex closure of the set of pure states {ip) such that there exists 
a pure state \£) G Sp with |(OV')l ^ f ■ The set Si j: is defined as 



Sij - {cmm I C e LOCCf , IV) e Spj}. 



(10) 
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It is not hard to see that the sets given in the Definitions are all compact, convex sets. In Appendix ^ we will 
prove that unlike the Schmidt number sets, these sets do not relate directly to positive linear maps. Nonetheless, the 
characterization of these sets in terms of a pure state and an LOCC^ operation is extremely simple. We are able to 
prove the following three theorems. The first is essentially a re-statement of a result by Jonathan and Plenio in 
the language of Si sets. 

Theorem 1 A density matrix p £ Si if and only if there exists a decomposition {pi, of p such that 

^Pi^i^i >- t^- (11) 

i 

Furthermore, the LOCC conversion 

^exact P, (12) 

is possible if and only if p £ Si . 

A straightforward corollary of this theorem is 
Corollary 1 Let p G Si , then E{p) < H{jl) where H{p) is the Shannon entropy of the vector jl and E is the 



entanglement of formation of p, 

E{p)^ min ^p,^;(|^,)(V'd), (13) 

where E{\^){ip\) is the von Neumann entropy of the reduced density matrix of \^) . 

The second theorem, dealing with the case of probabilistic conversions from pure to mixed states, combines the 
results in Ref. on exact conversions for mixed states with those of Vidal for probabilistic conversions between 
pure states 

Theorem 2 A density matrix p G 5*^^ if and only if there exists a decomposition {pi, IV'i)} of p such that 

pY^PrX^^)-"" pi. (14) 
i 

Furthermore, the LOCC conversion 

- P, (15) 
is possible with probability at least p if and only if p £ ^pp- 

For the last theorem, we need to introduce the notion of {p, /)-approximability. 
Definition 4 A density matrix p is (p, f)-approximable if and only if 

i 

for some decomposition p — \ipi){ipi\, and set of states {IV'i)} such that Tr^^ = Pi — (V-'ilV'i); ^^c^ 

E^i-^V- >- A*- (17) 



Remarks We could have associated a density matrix p' = G Si with the set of approximating states 

llV'i)}. It is not known whether our definition of {p, /)-approximability is equivalent with the more natural definition 
in terms of the square root of the transition probability F{p,p') = Tr(\/ p^/^p'pi/^); we could have called a state p 
{p, /)-approximable iff there exists a. p' £ Si such that F{p, p') > f. Because of the strong concavity property of the 
square-root-fidehty p| , 
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F{p,p')>Y,\{ij,m\ (18) 

i 

whenever p is (/i, /)-approximable according to Definition ^, it follows that p is also approximable according to a 
definition which uses the square-root-fidelity. Furthermore, even though there exists a pair of decompositions which 
saturate Eq. (|lj), see Ref. ||lj], these optimal decompositions need not correspond to ones which obey the rel ation 



of Eq. (y_7|). While we will further investigate the connection between F and (/2, /)-approximability in Section [ID, 
it remains unclear whether this second definition would allow for Theorem 3. 

The following theorem and its corollaries, when combined with the Vidal-Jonathan-Nielsen results [Q, create a 
framework for answering the question of approximating a mixed state using LOCC operations, starting from a given 
pure state. 

Theorem 3 A density matrix p £ Si j, if and only if p is {fl, f )-approximable. 

It turns out that for two qubits a simplification takes place: 

Theorem 4 For two qubits Si ^ ~ Sp^x where x refers to either exact, probabilistic or approximate conversions. Let 
S be the set of all bipartite 2-qubit density matrices. We have 

< g < 2^2 : %i,;,2),p=9 = (19) 



1 > g > 2^2 : 'S'(^i,^2),p=, = (20) 



Thus for two qubits, the probabilistic sets are identical to each other and identical to the exact sets. 
Theorems also give rise to necessary conditions for LOCC conversions from a mixed state to either a pure or 
mixed state. From Theorem |l| we can conclude that 

Corollary 2 The LOCC conversion 

Pi ^exact P2, (21) 

where pi is a mixed state is not possible if there exists a pure state \ipp) such that pi G Si , but p2 ^ Si . 

For probabilistic conversions, Theorem]^ leads immediately to 

Corollary 3 Suppose that pi G Si ^ and that, in particular, it is possible to obtain pi from \ipp) with some probability 
Pi > P- Now suppose p2 ^ ^p p- '^^^ maximum probability of success P2 of obtaining p2 from pi via LOCC operations 
satisfies p2 < P I Pi- 



ll. LOCC CONVERSIONS 



A. Preliminaries 

In proving Theorems d and I one observation turns out to be crucial. 

Observation 1 When a density matrix p can be obtained by LOCC from a pure state \i^p) with probability 
p G (0, 1] then there exists some decomposition {pi, of p into pure states, i.e. p = '^^Pi\^i){ipi\, such 
that we can obtain the ensemble of states j-i/'i) each with probability ppi by LOCC from the state \il)p). 

It is essential that the LOCC protocol starts with a pure state. Consider a unilateral operation by, say, Alice, on a 
pure state j-i/'): 

\i^)^{p..,sm{m, (22) 

where Si is some CP map, characterized by Kraus operators {^j.}- Instead of performing Si, she can perform the 
more 'fine-grained' measurement 
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(23) 

where = (g) and q,j = A^'' (g 1|V')(V'|, which gives her an ensemble of pure states. The extra bits of 
information that she obtains will be communicated to Bob in the next round. By replacing every 'mixing' operator 
in the LOCC protocol by a fine-grained measurement projecting on pure states, Alice and Bob obtain a new protocol 
which produces an ensemble of pure states as output, and the ensemble forms a decomposition of the mixed state 
that was obtained in the original LOCC protocol. 



B. Exact 

1. Proof of Theorem^ 

The proof of this Theorem follows quite straightforwardly from Observation 1 and the results of Jonathan and 
Plenio [0. Theorem 1 of Ref. states that we can obtain the ensemble {pi, IV'i)} (Si-Pi = 1) t>y LOCC from a 
state 14' ft) if and only if 

^p^X^, y A*. (24) 

i 

Now, if we can obtain the ensemble {pi, {ipi}} by LOCC, then we can obtain p = '^.iPi\ipi){''Pi\- On the other hand, 
if we can obtain p by LOCC from the state lipp), then we can obtain some decomposition of p by Observation 1, and 
thus for this ensemble Eq. (^) should hold. 

The equivalence between Eq. (^) for some ensemble {pi, \ of p and p G Si follows from inspecting the 
Jonathan-Plcnio protocol for converting \ijjp) to the ensemble {pi, It is a two step protocol: 

\lpft) ^exact 14) ^LOCC* P> (25) 

where \4) G Sp and its Schmidt vector is given by = X]j Pj'^V'j ■ Thus, every p for which there exists an 
ensemble obeying Eq. (p^), can be written as p = C{\4){4\) where C is some LOCC^ map and \4) € Sp. These 
observations prove Theorem pi. □ 



2. Consequences 

Note that since the majorization conditions treat Alice and Bob on an equal basis, we can alternatively write 

S+ = {C{\4){4\) I C e LOCC?, IV) e Sp}. (26) 

One might ask whether the extra complication incurred by defining the sets Si is necessary. While we will find later 

on that in the case of pairs of qubits the sets 5*^ are equal to the simpler sets Sp, the following example illustrates 
that for pairs of n-level systems, with n > 3, that the sets generally do not coincide. 

Example 1 Let \4a) ~ |00) and \4i2) = "^(IH) + |22)) then consider the density operator p = (1 — e)\4>o){4'o\ + 
e|V'i2)(Vi2|. ie<p= (l-e/2,e/2,0). Then for e e {0,1) , p e S+ but p <^ Sp. 

First, notice that A^^ — (1, 0, 0) while A^^j = (i, i, 0). Thus, (1 — e)A^|-, + eA^j^ = jl and we have that p e Si . 

Now we need to demonstrate that for any decomposition p ~ pi\(l)i){4>i \ there always exists an i such that A^^ )/- ft. 
Let 

Pi2 = (|l>(l| + |2)(2|)®(|l)(l| + |2)(2|), (27) 

and Pq — |00)(00|. Then P12P-P12 = e|Vi2)(V'i2| from which it follows that Pi2\4'i){4>i\Pi2 « |'0i2)(Vi2| and similarly 
Po\<t)i){(pi\Po oc |00)(00|. Together this implies that 

Ic/-.) = (Po + Pi2)|0.) = a.|00) + A|ii) + ^|22), (28) 

for some complex amplitudes ai,fii. Letting a*ai = ri and (3* fit — Si, we write A^^ = {ri, Si/2, Si/2), assuming 
ri > Si/2. If both r,; and Si are nonzero, then A,/,. )/- p. On the other hand, any decomposition of p into pure states 
each with either ri — or st — will necessarily contain terms for which \(j)i) = \4i2)- Thus, p ^ Sp. □ 
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C. Probabilistic 



1. Proof of Theorem^ 

First, we prove that if a state p has a decomposition p = j P j) {"^ j \ satisfying the condition 



P 



Y,Pj^^. (29) 



then p e p- Let us define the state with Schmidt coefficients A* = '^jPj^^.- Eq. ( p9|) directly imphes 

that £ Sp^p and thus by Vidal's result on probabilistic pure state conversions, Eq. we can obtain with 
probability p by LOCC from \ipp). Then an LOCC^ map as in the last step of Eq. (^) results in p. This shows that 
that p e si p (and also that we can obtain p by LOCC with probability p). 

On the other hand, let p G ^pp^ P = ^iWW) ^i^h C £ LOCC^ and 1-0) G This implies that we can 

obtain p with probability at least p by LOCC from li^p). By Theorem 1 we have that there exists a decomposition 
{pi, \ ipi)} of p such that 



and by Eq. (H), we get 



P 



'^Pi^ipi >- ^ipi (30) 



'^Pi^4>, > P^i! >^ P, (31) 



which is the desired condition. □ 

As before, this result allows us to write down an alternative characterization of the sets Si namely 

Sip = I ^ e LOCC?, 1^) e Sp^p}. (32) 



2. Relation with Schmidt number sets 

What is the relation of the sets 5*^^ with the Schmidt number sets For every vector p, there is some probability 
p by which \ipp) can be mapped onto \'^k) where k is the Schmidt rank of p. Therefore Sk C Si^ where pmax is 
given by the maximum probability p such that pXxj,^ y'^ p. Note also that Si^^ C Si ^ where q < Pmax- On the 
other hand, we have Si^ C Sk for all p > 0. This follows from the fact that (1) jV'/j) has Schmidt rank k and (2) any 
state \ipf_i) can be reached by LOCC from jvl^fe) with probabihty 1, due to the condition in Eq. (^), and (3) from \ipp) 
all states in can be obtained by LOCC with probability at least p. 

Therefore 

Sk ^ Si^^, yO < q < Pmax, (33) 

that is, for small probability up to some maximum, the sets Si ^ collapse onto the set 5*^ and onto each other. 



D. Approximate conversions 

The following definition will be useful in what follows. 
Definition 5 (Optimal /i-approximation) A density matrix p' — J^i Wilii^'ilr where Pi = and 

"^Pt^ii', >- (34) 

i 

is an optimal ^-approximation to p — \tpi){'ipi\ when 

i 

is maximal with respect to all decompositions of p and sets of states jV'i) for which condition Eq. holds. 
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1. Proof of Theorem^ 



Let p be a (/I, /)-approximable state and p' be an optimal /i-approximation to p. Thus there exist decompositions 
P^E^ and p' - IV'^X^^I such that 

Y.\(^^\^'^\^f- (36) 

Since p' G S^, we can find unitary operators U'^ and V/, non- negative diagonal Kraus operators Ai = 
diag(^Ai4, . . . , y^Xi^n) satisfying J2i^f = 1 and a state |V'p') e % such that = (f// (g) y/^i)|?/'p')- Simi- 
larly, we can find unitary Ui and Vi, operators Bi = diag(^/ii_i, . . . , y/Jti^), X^i — ^ and a state ji/jp) such that 
l^/'j) = {Ui ViBi)\tpp). Furthermore, the Plenio- Jonathan protocol requires that jV'p) and \ipp') be in Schmidt form 
in the bases that diagonalize their respective Kraus operators Ai and Bi, which in both cases is the standard basis. 
Therefore, IV-p') = J2k ^) I^p) = Sfc VPk\k, k). 

By Lemma 1 of Ref. M, the absolute value of the inner product between two states 1-0) and with given Schmidt 
coefficients is maximized when the states are chosen to have the same Schmidt basis. We can therefore assume that 
U'^ = Ui and V/ — Vi since p' is an optimal p-approximation to p. We can then calculate 

E I 1^^) I = E ^o^kKkVP^k- (37) 



However, since ~ "l^it^i-k ~ 1; know that y\j4M^ < 1- Therefore, again using the optimality of p', 

we can conclude that pi^k = K,k- It follows that 

m^p')\-J2\(^^\^'^\^f- (38) 

i 

Therefore \ipp) g S'^,/- We have shown that when p is (p, /)-approximable, we can write p — £(|'0p)(V'pl) where 
\ipp) £ Spj and £ is an LOCC^ map. 

On the other hand, let p £ Si j:] we can write p = C{\^p){ip\) where lip) £ Spj. Consider an optimal p-approximation 
IV'') £ Sjx to \^p). As an approximation to p we take p' = £{\^p'){ip'\). We write p = J2i \''Pi){^i\ and p' = J2i \'4^i){'^i\ 
where {Uii^ Ai)\'ip) = [ipi) and similarly for {ipi), where Ai are the Kraus operators of a local superoperator implemented 
by Alice. Then we have 

E 1(^*1^^)1 = E K^I^A. ® > KV-IV-')! > /, (39) 

i i 

where we used X- |a;i| > | J2i Therefore p is (p, /)-approximable. □ 



2. Connection to square-root fidelity 

We will now investigate in more detail the relationship between square-root fidelity and (p, /)-approximability. For 
convenience, we begin with 

Definition 6 (p-fidelity) Let p and a be density operators. The jl-fidelity is defined to be fp{p, a) = max^. 

where p = J2i \4'i){(t>i\ o.'^'d a — J2i with J^iPi^i'i A* for Pi = (ipil'ipi)- If no such decomposition exists then 

set ffi{p,a) = 0. 

Furthermore, define the following two optimal fidelities: 

Frr,ax{p, p) max{F(p, a):a£ Si} (40) 
fniaxip, p) = max{/^(p, a) -.(T £ si}. (41) 

In this section, we will develop some properties of Fmax and fmax that might shed some light on the relationship 
between them. We can concisely state our results on approximation by noting that Fmax (p, p) > fmax (p, p) and 
p £ si j: if and only if fmax{p, p) > /• However, since Fmax is a more natural definition of the optimal fidelity of 
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approximation, we would ideally like to re-phrase the second result in terms of Fmax- We will not be able to do so 
but can offer some suggestive partial results. 

To begin, both functions are easily seen to be consistent with the partial ordering on Schmidt vectors: X fl 
implies that Fmaxip, X) < Fmaxip, fl), and similarly for fmax- Furthermore, notice that for a pure state 10) and 
Schmidt vector fl, 

frnaxiimip) = FmaxUmip)- (42) 

One way to see this is to recall the result from Ref. that the best approximation with respect to F to a given pure 
state by a mixed state a G Si can always itself be taken to be given by a pure state lip) e Si . The result then 
follows because fp{\(l)), IV')) and F{\(t)), \tp)) agree if A|^) >- fl. 

The functions are similar in other ways as well, both them obeying a strong form of joint concavity. Let us begin with 
Fmax- Suppose that Fmax{Pj,Pj) = F{Pj,o-j) for some aj G Si,. If we set /2 J2j QjPjy then we find Y^j Qj'^j ^ Si 
so that 

FmaxC^PjPj,^ qkPk) > PC^PjP3,Yl 'i^^^^ (^^3) 

3 k j k 

>Y.Vm3F{P3^^3) (44) 

3 

= 51 VmjPrnaxipj.P])- (45) 
3 

Let us now move on to the joint concavity of fmax, for which the proof is similar. Suppose that fmax{pj , Pj) = 
fpjiPj,<Jj), where aj € Si,. Likewise, assume that the optimal decomposition is given by fpj{pj,crj) = J2i \ {4>ij\''Pij)\, 

where J2^'^i3^i'^j >~ flj f^'' = (V'ijlV'ii)- Then, since Qj'^is^ipij >~ J^j Ist^j' ^^'^ 

fmaxiY^ Pj Pj ,^qkPk)>Y^ VP3^\ ('I^v I V'y- )\=Y1 -^aZ, (Pj ' ) = XI Z™"^ (^J ' A'j ) ' (46) 

j k ij j j 

which is the desired inequality. 

Summarizing then, for a function g equal to Fmax or fmax, we have 

1. (Consistency with partial order) X)^ fl implies g{p,X) < g{p,p). 

2. (Agreement on pure states) For a pure state \4>) and Schmidt vector /I, g{\<P'){4>\, p) is equal to the maximum 
over \^) &SpofM^)\. 

3. (Joint concavity) If p ^ J2iPi\(t>i){'Pi\ and p J^j I3P3 then g{p, p) > J2i VPiOidiPi, Pi)- 

There is a difference, however. For fmax we also know that given p and p, we can find a decomposition p = 
Z]iPi|0i)(<^j| and an ensemble {g,, with J2i >- P such that fmax{p, p) = ^/PiQi\{<l>i\'^i)\- Since \ {(t>i\iki) \ < 
frnax{\4'i){4'i\,^ipi) by definition, summing inequalities gives 

fmax{P,P) < X VPiOifniax {pt,X^,). (47) 
i 

Combining with joint concavity then implies that 

fmax{p,P) Xx/Pig7/ma2;(|0i)((/'i|,A^J- (48) 

i 

Thus, fmax is the smallest function consistent with the three properties listed above. 

This is quite a curious situation. The results of Ref. [ p^ immediately imply that the square-root fidelity F{p, a) 
(note that the arguments of F are different than those of Fmax) is the smallest function g satisfying the following 
different list of properties: 

1. (Agreement on pure states) For normalized |(/)) and g{\(t)){(j)\, |'0)(V'l) = l(0l'0)l- 

2. (Joint concavity) li p = J^iP^Pi and a = J^j ^3(^3 then g{p, a) > J2i \/PtQi9{Pi, o'i)- 

Thus F arises as the smallest function satisfying a set of conditions which appear very similar to those conditions 
satisfied by Fmax and fmax- It seems plausible, therefore, that Fmax, which is defined in terms of F^ might itself be 
the smallest function consistent with properties 1-3 of Fmax and fmax- If that is the case, then fmax — Fmax- 



8 



3. Applications 



Theorem 3 does not tell us how to determine the maximum /i-fidehty with which wc can approximate a density 
matrix p, given a starting state To answer this question will generally require a further optimization over 

decompositions of p. Let us consider an example. 

Example 2 Let = be a maximally entangled bipartite state of Schmidt rank m and let 

Am be the Schmidt vector of |4'm)- The best approximation to a state p achievable starting from using LOCC 

operations will have an optimal approximation with 

fmaAp, A,n) = (1 - (49) 

where Ei{p) is one of Vidal's entanglement monotones J^, defined as 

Flip) = nnnY,P^Eli\A){^^\), (50) 

i 

in which £'((|'0i)(V'i|) = X]fe=i A^. and £ = {pi, {tpi)} is an ensemble realizing p. 

In Ref. 0, it is demonstrated that if |0) = J2^=i\^\''-A)\iB), with f3i > Pi+i then the best approximation to |0) 

achievable using LOCC operations and starting from the state \'^m) has square-root-fidelity fmax — (X^I^li A)^^^- 
Theorem 3 shows that the optimal approximation for mixed states is that derived from some Jonathan-Plenio precursor 



state of the form y X^jPiA^Jiyi)!*^), where p = J2j Pj\^j)i'^j\- Therefore, the best approximation to the mixed 
state p has fidelity 

fmaxip, Am) = max Y^pj ^ A^ I = min^Pj ^ A^,^. , (51) 

\ j j=l / \ j i=m+l J 

where the optimizations are to be performed over all decompositions of p. The last line of the above equation, however, 
is simply (1 - £;m+i(p))i/2. □ 



E. Two Qubits: Proof of Theorem g 

1. Exact 



First we relate the set Sp to the entanglement of formation: 
Proposition 1 The entanglement of formation E of a density matrix in C2 ® C2 is 

E{p)= min H2{p), (52) 

where H2{.) is the binary entropy function, i.e. H2{x) — —a; log a; — (1 — a;) log(l — x). 

Proof The Wootters formula for the entanglement of formation [|l3| tells us that in the optimal decomposition 
{pi,\ipi)} of p, every state Itpi) has an equal amount of entanglement. Therefore if p has entanglement E(p), the 
density matrix p must be contained in the set Sp where H2{p) = E{p). Conversely, ii p € Sp then E(p) < H2{fl).0 

From this Proposition it follows that Sp must be closed under LOCC and therefore Sp = Si . Furthermore, for two 
qubits we know that the partial order on the pure states induced by majorization is a total order, which is characterized 
by a single parameter, the smallest eigenvalue fj,2 of the reduced density matrix of the pure state. Therefore we can 
characterize the vector p with this single parameter p2 and it follows that ■^'^^ when p2 < P2- 

As a corollary of this result we obtain that we can convert jV'/i) to p by LOCC if and only if p e Sp, which means 
that E{p) < H{fl). Therefore the minimal entanglement costs for preparing a single copy of p exactly is E{p). This 
result has been independently found by Vidal |15|] . 
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2. Probabilistic 



Let us now consider the probabihstic sets and prove that for 2 qubits Si^ = Sjx^p. First we note that the relation 
with the Schmidt number sets gives us the following: 

S = 82 = 5'(^i,/:i2),P=9' (53) 

when q < Pmax = 2/U2- Here ^2 is the Schmidt number 2 set, which is identical to the set of all bipartite two-qubit 
density matrices, S. Now let q > 2^2- A density matrix p G Si ^ if and only if there exists a decomposition {pi, \tpi)} 
of p such that 

Q^^PAli < 1^2, (54) 

i 

where A^. is the smallest Schmidt coefficients of \ipi). Since P2/Q < 1/2, this condition is identical to the requirement 
that p € 8^'=(n2/q i-n2/q)' ~ % ^'^^ qubits, this condition is again identical to requirement that there 

exists a decomposition {pi,\il>i)} of p such that for all i 

9A^, < M2, (55) 

which implies that p G S^^g. So we have shown that for q > 2p,2, the exact sets 5'(/i2/g,i-/*2/9) ^^'^ probabilistic 
sets S(^fjt^^^^)^q are, in fact, identical. 

3. Approximate 

Finally, let us move on to the approximation sets. We'll begin by considering pure states. Suppose that ^ Sp, 
and let S Sp such that |('0|'0')l = fmax is optimal. The condition that \ip) G Sp^j is that fmax > /• Suppose now 
that IV') has Schmidt coefficients (a, 1 — a), where a> 1/2. It is clear that fmax is a strictly increasing function of a 
on the interval [1/2, jii], from which it follows that fmax > f, or, equivalently, that jV') e if and only if a > ^[ 

for some p'l £ [1/2, ;Ui]. If we set p'2 = 1 — fi'i, it follows immediately that Spj' = Sp' since the extreme points of 
these convex sets are pure states. It then follows from the definitions that Si = Si,. We saw above, however, that 
S+=Sp,. Thus, S+j = Spj. 

III. CONCLUSION 

We have provided a unifying framework for exact, probabilistic, and approximate LOCC conversions from pure 
states to mixed states. In each case we have found criteria defining exactly when these conversions are possible, with 
the caveat that the criteria are always expressed in terms of the existence of 'optimal' (lec;ompositions of the target 
mixed state having some easily verified property. This work does not address the question of how to find these optimal 
decompositions of a density matrix p for exact, probabilistic or approximate conversions of an initial state jV'/i) into 
the state p. This problem could be as hard a determining optimal decompositions for the entanglement of formation, 
but may well be simpler. In addition to resolving that question, with our framework established, a host of open 
questions present themselves. The most pressing is the question of determining whether the function fmax{p,fj) we 
have been studying here is or is not equal to the more natural Fmax if, p) ■ Another question of interest is the relation 
between the sets Si and the optimal decomposition of p with respect to its entanglement of formation. 
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APPENDIX A: POSITIVE LINEAR MAPS 



We may ask whether, say, the sets Sp or St are characterized by positive hnear maps. Let us define a /i-positive 
hnear map in the foUowing way, as a generaUzation of fc-positivity. Here /2 is a n-dimensional Schmidt vector. We 
would propose the foUowing definition. The Hnear Hermiticity-preserving map C is positive with respect to {ipp), or 
/i-positive if [1 (E) C) (p) > for aU density matrices p € Sp. It turns out this impHcs that C is fc-positive where k is 
the Schmidt rank of the vector jl. The property of /2-positivity imphcs that for aU |?/;) G Sp and for arbitrary |0) we 
have 

(0|(i®/:)(|V')(V'l)l0> >o. (Ai) 

We can always write = (A (g) 1)|$) where |$) is some maximally entangled state. By commuting A through, we 
obtain 

mi^c){\^A){^A\m>o, (A2) 

where Ii/j^i) = {A^ (g) There is always an operator A and therefore a state and |$) such that is 

proportional to the m axi mally entangled state with Schmidt rank fc, where k is the Schmidt rank of {tp) itself. Then 
the condition in Eq. ( [A^ ) becomes the condition for fc-positivity of C as given in Ref. 1^. Thus jl- positivity implies 
fc-positivity, whereas Sp C Sk- 
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